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Motivated by the problem of domino tilings of the Aztec diamond, a 
weighted particle system is defined on N lines, with line j containing 
j particles. The particles are restricted to lattice points from to N, 
and particles on successive lines are subject to an interlacing constraint. 
It is shown that marginal distributions for this particle system can be 
computed exactly. This in turn is used to give unified derivations of a 
number of fundamental properties of the tiling problem, for example the 
evaluation of the number of distinct configurations and the relation to 
the GUE minor process. An interlaced particle system associated with 
the domino tiling of a certain half Aztec diamond is similarly defined and 
analyzed. 

1 Introduction 

The analysis of certain tiling models is of common interest to both combinatorics and 
statistical mechanics. As an explicit example, consider a so-called (a, b, c) hexagon. 
This is a hexagon with integer side lengths a, b, c, a,b,c — side lengths a vertical by 
convention — reading anti-clockwise, and all internal angles ^ (see Figure [I] for an 
example). Such a hexagon can be tiled using three species of rhombi, each with side 
lengths 1 and angles ^, 4^. The three species of rhombi are distinguished by their 
orientation - down sloping, up sloping or neutral in slope, reading left to right. As 
illustrated in Figure [l] it is immediately clear that a particular tiling of the hexagon 
can be uniquely specified by a family of non-intersecting lattice paths. These all start 
and finish one unit apart, and move up or down half a unit at each step (reading left 
to right). 

In this article our interest is in tilings of the Aztec diamond by 2 x 1 dominoes. The 
Aztec diamond of order N is the union of all lattice squares within the diamond shaped 
region {(x,y) : \x\ + \y\ < N + 1}, and the dominoes may cover the lattice squares 
by being placed horizontally or vertically. As with the hexagon tiling of the previous 
paragraph, such a tiling of the Aztec diamond can be uniquely specified by a family of 
lattice paths. To see this, with the top left lattice square specified as white, introduce 
a checkerboard colouring of all the lattice squares making up the Aztec diamond. For 
a horizontal domino which covers a white-black (black-white) pair of squares when 
reading left to right, no segment (a horizontal segment) of path is marked. For a 
vertical domino which covers a white-black (black-white) pair of squares when reading 
top to bottom, a right-up (right-down) segment of path is marked. This results in a 
family of N non-intersecting lattice paths, with segments up sloping, down sloping or 
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Figure 2: The 5 non- intersecting paths corresponding to a particular domino 
tiling of the Aztec diamond of order 5. 
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Figure 3: An example of the shading of an Aztec diamond of order 10, rotated 
45°. Here, the lines pass through the black squares of the checkerboard colour- 
ing, and we can see that all the E and S type dominoes (the shaded dominoes) 
are intersected by lines in their left half, providing an easy way to check the 
shading. 

horizontal, starting at equally spaced points on the bottom down sloping edge, and 
finishing at the corresponding points on the bottom up sloping edge. See Figure 2 for 
an example. 

In the case of the tiling of the hexagon, Figure[l]makes it clear that complementary 
to the non-intersecting paths, the tiling configuration can equally as well be specified 
by recording only the centres of the neutral in slope rhombi. These centres can in turn 
be regarded as a particle system on the set of vertical lines naturally associated with 
the (a, 6, c) hexagon. This particle system has the peculiar property that the number 
of particles equals the number of lines for lines 1,2, ... ,b, then stays constant for lines 
b+ 1, . . . , c, and then equals b— 1, 6 — 2, . . . , 1 for lines c+ 1, . . . , c+6— 1. Furthermore, 
the particles are constrained by interlacing constraints, the details of which are evident 
by inspection of Figure [T] 

Our interest is in the interlaced particle system implied by a domino tiling of the 
Aztec diamond. For its specification, with the Aztec diamond checkerboard coloured 
as already described, let the horizontal dominoes such that the left square is colour 
black (white) be called of E (W) type. Similarly, let the vertical dominoes such that 
the top square covered is black (white) be called of S (N) type [6]. Suppose now 
that the E and S type dominoes are shaded and numbered lines added (see Figure 
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Figure 4: An example of two different tilings with the same particle picture. 
As can be seen, where a particle is not adjacent to another particle on the next 
line, a 'square' is formed, that can be tiled in two different ways. 

3) . Each line k passes through the interior of k shaded tiles, and these intersections 
are considered as specifying the positions of k particles [151 116| . In an appropriate co- 

(k) (k) 

ordinate system, these particles occupy distinct positions x\ > • • • > x k restricted 

to the lattice points 0, 1, 2, . . . , N on line k (k = 1, . . . , N). Most importantly, the 
particles must satisfy the interlacing condition 

4% < < 4 k) for i = l,...,fe-l (1) 

A crucial point in relation to our study is the inverse of this mapping. Consider 
co-ordinate a;' fc_1 ' on line k — 1. Suppose furthermore that the interlacing condition 

holds with strict inequalities. Then there are precisely two domino orientations 
corresponding to . On the other hand, if either inequality in |l]) is an equality, 

there is just a single possible domino orientation corresponding to -1 ' (see Figure 

4) . Importantly, this means that unlike with the hexagon, given a random tiling of 
the Aztec diamond with every possibly tiling equally likely, the corresponding particle 
system must be weighted. 

We remark that beyond the theory of tilings of the hexagon and the Aztec diamond, 
interlaced particle systems with varying numbers of particles occur naturally as the 
eigenvalues of successive minors of Hermitian random matrix ensembles [T] [5] 111)] . 
In fact we will show that in a certain scaling limit the particle system for the Aztec 
diamond converges in distribution to the eigenvalue process for the minors of Gaussian 
complex Hermitian random matrices (minors of the GUE ensemble). Also, although 
not a theme addressed here, we remark that interlaced particle systems are a rich 
source of determinantal point processes [211 [TT1 l22l [T6l l3l 131 l4l l2l ITSl \W\ . 

In this paper, we use the underlying particle system to rederive fundamental results 
about random domino tilings of Aztec diamonds, including the number of possible 
tilings, the multi and single line probability density functions (PDFs) for the positions 
of shaded particles, the limiting large N shape of the disordered region (Arctic circle 
effect), and the relation to the GUE minor process from random matrix theory. We 
consider a different particle system corresponding to the domino tiling of a half Aztec 
diamond, and exhibit analogous properties, in particular the limiting large N shape of 
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the disordered region (now half the Arctic circle) and the relation to the anti-symmetric 
GUE minor process. 



2 One and multi-line PDFs 

Consider a sequence of vertical lines in R 2 , with the fc-th line at x — k and containing k 
particles. Let p(x ( - m \ . . . , x^ n ') be defined as the joint probability that the i-th largest 
particle on line k is at (k, x^) for k = m, . . . ,n and i — 1, . . . , k. For the particle 
system relating to a random tiling of an Aztec diamond of order N, we know from the 
discussion about jl]) that x\ h ^ must obey the restrictions 

< < xf < N, 



x 



(fc+i) J*) ^ ^(fc+i) 
i+i 



< x\ K > < x\ K+1 > , (2) 



x\ h) e z, 



(the second of these is just |TJ with k i-t k+1). We also know that although each tiling 
is equally likely, each particle system is not, since a particle system is not uniquely 
defined by a tiling. To account for this, we introduce the notion of adjacency. Let a;^ 



be called adjacent if for some j, x^ = Xj +1 \ Each particle that is not adjacent (and 
is not on the last line) represents a tile with two possible orientations, and so must 
be weighted by 2. Given that there are ^N(N — 1) particles not on the last line, the 
joint PDF for the entire particle system is given by 

r,JV(JV-l)/2-a(:i;< 1 ),..., : E (JV - 1 ») 

p(x W , ■ ■ ■ ,* (JV) ) = - A X(* (1 \ ■ • • ,* (A °) (3) 

An 

where a, the number of adjacent particles, is given by 



a{x {m \...,x {n) ) = n{x (k) ) = 5 x f\ x f+V + S x f\ x f+V ( 4 ) 

k—m i=X,...jk 
k — m,...,n 

An is the number of possible tilings of an Aztec diamond of order N, and 

f (m) <Wh / 1, if {x {m \...,x {n) ) obey pi 

^ ) = | 0, otherwise U (5) 

Proposition 2.1. For the particle system corresponding to uniform random tilings of 
the Aztec diamond of order N 

I fm) <N)\ A(x^ m ') , ( m ) o\n. . . 



where 



and 



a(^)= n (^-^f) m 

l<i<j <n 



m — 1 



^ = a JV 2- jv ( jv - i »/ 2 n « (8) 

Proof. The m = 1 case is true from pi. Assume the m = n case is true. Then 

D (a; (n) x iN) ) - Afz'"') y (x (n) x {N) ) (91 

P[X '■•■' a; ' DnN2a (^),...M N - 1 '>) X( ' (> 
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Summing on the n-th line gives 



p(x 



„(™ + l) 



E 



E 



1 ~~ 2 

x(* ( " +1) , 



A(a: ( " ) )x(a: ( " +1) ,...,^ 

D n , JV 2-(- < " ) >-.- (lv - 1) ) 



-Dn,jv2 a (- ( " +1) .--. a:(JV " 1) ) 



det 



V — 



i,j — 1 , . . . ,n 

(10) 



where we have set at — 



„(«+!) 



„(«+!) 



The sum in the determinant is a 
polynomial lunction of a% and bi with highest degree term (6| — o|)/ j. Since the lower 
degree terms will have the same dependence on ai,bi for each row i, they can be 
cancelled out by column operations. Thus 



p(x 



(n+l) 



J D njJV 2"(- ( " +1) .-.- (JV - 1) ) 
A(x (n+1) ) 



det 



(n + l) 



i,j — l,...,n 

(JV) 



,^ VJ ) (11) 



n!A l ,iv2-(- < " +1) .-- <JV ™ 1) )' 

where the determinant evaluation follows by noting that it must contain A(x' n+1 ') as 
a factor, and is of the same degree as A(a;' n+1 '). The case m = n + 1 has thus been 
established, provided Dn+i.jv = n\D n ,N, which is indeed a property of IpJ. □ 

To find An, we introduce virtual particles {xj }i=i,...iv+l to the system with 
the requirement that also obey Q. We note that the only possibility is x^ N+1 ' = 
{0, 1, . . . , TV}. Following the method of the proof of Proposition 2.1 but beginning 
with the PDF 



p(x 



(i) 



- ) N(N + l)/2-a(x( 1 \...,x {N '>) 



At 



(i) 



JJv+i)n 



(12) 



(the term N(N + l)/2 in the exponent results from there now being TV + 1 lines), we 
end up with the single line PDF 



p{x(N+1)) = 



But x_ 

p(x^ = {0,l,...,N}) 
Hence we conclude 



N + 1 — i, which substituted in |13j) gives 

9 JV(JV+l)/2 TT 
* lLl)<i<j<N 



(13) 



(14) 



A N = 2 N(N+1)/2 (15) 

The result ( 15 1 for the number of domino tilings of the Aztec diamond was first derived 
by [6]. Since then a number of derivations distinct from those given in [6j have been 
found, for example [141 [7] • The present derivation using the particle picture appears 
to be new. 

It remains to compute the single line PDF, which is gotten from (JsJ) by summing 
over x (m+1 \...,x (N) . To perform the summations we introduce a second set of par- 
ticles representing the unshaded tiles. As with the shaded tiles, we want the fe-th 
line to have k particles, so we label the lines from right to left in the y picture. Because 
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every position must have a shaded or unshaded tile, and no position can have both, 
the j/*' are defined such that 

x (n) [J y (N+l-n) = {0,1,... N} 

x (n) n y (N+l-n) = 0> (16) 

Because the unshaded tiles, when viewed right to left, obey the same probabilis- 
tic law as the black tiles view left to right, the formulas for p{x i - Tn \ . . . , x^) and 
p(t/ m ' , • • • , j/ n ' ) are the same for all m,n. 

We now express (J©J) as a function of the yf\ To begin, using the fact that 

N 

A({o,i,...;v})=n i! 

i=l 

we have 

It remains to calculate 2 a ^ ''' " ,:E< ^ in terms of y\ k \ 

Proposition 2.2. Consider two lines n, n+1 in a particle system as defined above, but 
generalised so that there are N* possible positions for particles on each line. Let this 
system be filled with x and y particles, such that every possible position has either an 
x or y particle, and no position has both. Furthermore let the lines labels be changed 
to n* and n* — 1 respectively when considering the y particles. If line n has a x- 
particles {x^, . . . , a;!"'} (and therefore N* — a y-particles {y' n \ . . . ,y^,^_ a }) and 
line n + 1 has b x-particles {x^ l+1 ^ , . . . , x^ +1 ^} (and therefore N* — b y-particles 
{y[ n _ , . . . , yjy.Zfc } )then, for a defined as above, 

a{x (n) ) = a{y {n *~ 1) ) + a + b-N'" (19) 

Proof. There are a x's on line n. Of these, a(x^ n ') are adjacent. Therefore, exactly 
a — a(x^) of the :r's on line n are not adjacent. Noting that line n + 1 is the lower 
numbered line in the y picture, this means that exactly a — a (x {n) ) of the j/'s on 
line n + 1 are not adjacent. Since there are, by definition of a, N* — b — a(y' n _1 - ) ) 
non-adjacent y's on line n+1, 

a - a(x (n) ) = N* -b-a{y {n *- iy ) (20) 

and so (Tl9l) follows. □ 



Using Proposition |2.2| with a = n, b = n + 1, N* — N + 1 and n* = N+1 — n gives 
a{x (n) ) = a(y (N ' n) ) + 2n- N (21) 
a(x (m \...,x (N - iy ) = a(y W ,...,y {N - m) ) + (N -m)(m-l) (22) 



So, applying ( 18 1 and ( |22[ ) to (|6| we have 

p(v W „(nK niTo 1 ^-*)! A(^)) x (yW,...,^) 

But we knowp(a; (1) ,...,i w ) = p(y m , . . . , y {n) ). Finally, using the same inductive 
method from the proof of Proposition |2.1| we end up with 

This has been derived using different arguments in [13] (in particular the weighted 
particle system is not specified by ) , where it is recognised as a particular example 
of a discrete orthogonal polynomial unitary ensemble based on a the Krawtchouk 
weight with p = 1/2) (see the Appendix). 
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3 Large N limits 



In 16 the weighted particle process corresponding to an Aztec diamond tiling, defined 
through its correlations and restricted to the first n lines, was shown in a certain scaling 
limit to coincide with the minor process for a certain ensemble of random matrices. 
These random matrices are the ensemble of complex Gaussian matrices Xwith measure 
proportional to e~ TlX ^ 2 , to be denoted GUE* (conventionally the GUE ensemble has 
measure proportional to e~ TrX ). The minor process is formed out of the correlated 
eigenvalues U" =1 {z ( --''}, z' 3 ' = (zi, Z2, ■ ■ ■ , z n ) denoting the eigenvalues of the j-th 
minor. This is known Q] to have joint PDF 

1 n n — 1 

l ne -(W/ 2 Yl ( 2 W-4"))H^ +1) >^) (25) 

" 1=1 l<j<fc<n 3=1 

where, with \a the indicator function for the set A, 



x{z u + D >z ^ ):= ^ 



and the normalization C„ is given by C n = (27r) n/2 . We can show directly that that 
the joint PDF for the weighted particle process tends to |25l in an appropriate limit. 



Proposition 3.1. Let the points Zf := (2y^ — N)/y/N be a rescaling of the points 

M 
n ' 

yf' have PDF p as described in (23 1, one has 



y\ , where N is the order of the Aztec diamond as described above. Given that the 



p(y (1 \...,y (n) ) -^p*{z {1 \...,z {n) ) as N -> oo (26) 



where p* is the PDF for the GUE* minor process as specified by (25 I. 
Proof. Let y = (z^N + N)/2 := g{z). Then 

j n 

p(y«, . . . = p ( 5 ( 2 (1) ), . . . , </(*<">)) J] II 9'(4 d) ) (27) 

i=l 2=1 

Clearly, 3'(^ (j) ) = y/N/2, so 

n s'(^ 5 ) = (f ) (28) 
<=i,...j 

3=1,. ..,n 



We now wish to compute p (g(z^), . . . , g(z^)j in the limit TV — > oo. Applying forms 
of Stirling's approximation (for large TV) 

(a/V + 6)! ~ %/2^^(a/V) aJV+ V aJV 
(aN + bVN + c)\ ~ %/2^aiV(aiV) aJV+,,y3?+c e i,2/2a - aJV 



to (|23|, and using (|27|) we have 

Z 

^27r)"/ 



A / (n) \ n — 1 

p(y« . . .,„<»>) ~ ^7^- ^ )2 II > (29) 



3 = 1 



which is (251. □ 
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We would also like to compute the region of support for this particle system. In 
the Aztec diamond tiling this corresponds to the boundary of the disordered region. 
The region of support on any given line n is the interval [a„,6„] in which, for a large 
enough number of particles, all the particles will lie within that region with probability 
1. Since we are taking n, the number of particles, to be large and n < N we must 
take N, the number of lines, to be large also. It thus makes sense to scale n so that 
the 'line label' s = n/N is a real number in [0, 1]. The region of support of the system 
will be the areas in between the graphs of a(s), b(s). To calculate the functional form 



of these boundaries from (24 1, one approach would be to use the fact that this PDF 
relates to the Krawtchouk ensemble. The necessary details have been given in [12] . 
Here we give a more physically motivated derivation, based on a log-gas picture [8]. 
The Boltzmann factor for a log-gas of N p particles has the form 



ii i**-*/]!' 

l<i<j<N p k=l 



(30) 



where /3 denotes the inverse temperature and V(x) is a one body potential, due to 
background charge density —p(x). Explicitly, 



V(x) := / p(t) log \t - x\dt. 



(31) 



A hypothesis of the log-gas picture is that for large N p and to leading order the particle 
charge density and background charge density cancel, so that the particle density is 
to leading order equal to p(x). 

In the cases that p(x) is supported on a single interval [a, b] (which we expect for 
the log-gas interpretation of ( 24 1 ) , normalization of the density requires 



p(t)dt = N p . 



(32) 



Furthermore, the explicit form of p(x) obtained by solving the integral equation (31 1 
is known in terms of V(x), and the boundary of the support is determined by the 
equations [8] 



V'ji) 

a y/{b-t){t-a) 

b tV'(i) 
a ^{b - t)(t - a) ' 



dt = 



dt = TlNr, 



(33) 
(34) 



As written, ( 24 \ is a lattice gas variant of the log-gas ( 30 1 in the case /3 = 2. In the 
limit n — > oo, the lattice gas approaches the continuum log-gas upon the substitution 



„(«) 



= Nt 



(n) 



(35) 



where, to leading order in N, < < 1. In terms of the co-ordinate t\ n> = ti, the 



(n) 



one body factor in ( 30 1 reads 



-2V(t) 



{Ni)\{N -Nt)\ 



(36) 



and recalling s = N/n shows N p = Ns. Thus solving (33 1 and (34 1 in the limit 
N — > oo gives (a(s),b(s)), the support in the variable t. 

From Q we know that a,b £ [0, 1]. Noting that V(t) = V(l — t) and inserting this 
into ( |33[ ), we have that for some c 6 [0, a = | — c, b = | + c. Computing 



lim 

N— >oo 



N 



log 



(37) 
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leaves us with 

,fc ilog(t/(l-t)) 



= dt = 2tts (38) 



= \l{\+c-t)(t + c-\) 
to solve for c = c(s). The change of variables t = | +lt leads us to a more managable 

fM±M* = M . (39) 

The integral can be computed exactly (a computer algebra package was used), giving 



1 - sfl - 4c 2 = 2s. (40) 

Recalling that c £ [0, |], we see that this has a solution for s 6 [0, |] only. For 
s £ (!)!]> physical interpretation of the relationship between s and c (increasing s, 
the number of particles, must not decrease 2c, the size of their support) leads us to 
define c(s) := | for s G (§,!]• So we have 



a(S) -\0 a e(|,l] 



(41) 



b(s) = {p+v^n^iF) . | [0,11 (42) 

Noteworthy here is that, because of the mirrored nature of the shaded and un- 
shaded tiles, the area of support of the unshaded tiles is related to the area of support 
of the shaded tiles by 

a a haded(s) = QWshaded(l ~ s), &ahaded(s) = &unshadcd(l — s), (43) 

so the disordered region of the Aztec diamond tiling, the area that has both shaded 
and unshaded tiles, is a perfect circle — the Arctic circle — in the limit N — > oo |12j . 



4 The half Aztec diamond 

Consider an Aztec diamond of order N = 2(M + 1) rotated by forty five degrees as 
in Figure [I] Define a restriction on the tiling of this Aztec diamond such that in the 
particle picture as defined above, a particle at x on line j implies no particle at x on 
line N + 1 — j. Because of the interlacing restriction, this means that in the tiling 
picture the whole middle column between lines k = M + 1 and M + 2 will consist of 
squares formed from a pair of dominoes rotated 45°. If we delete all these squares 
we are left with two halves. We will call these half Aztec diamonds of order M. The 
present half Aztec diamond model bears some resemblance to the Aztec diamond with 
barriers introduced in |20| . 

By construction, the tiling corresponding to two half Aztec diamonds are mirror 
images. We will call any tiling of an Aztec diamond of order N = 2(M + 1) formed 
from two half Aztec diamonds symmetric. With Cj^ the number of symmetric tilings 
of an Aztec diamond of order N and Hm the number of tilings of a half Aztec diamond 
of order M, we therefore have 

CV+i) = H M 2 M+1 (44) 

Here the factor of 2 M+1 corresponds to the number of tilings of the deleted squares. 

We would like to use the particle picture to compute Hm- We begin by noting 
that the joint PDF for the weighted particle system is 

r ) M(M + l)/2-a(x { - 1 '> ,...,x (M) ) 

P(x (1 \- ■ • ,- (M+1) ) = 2 = ^ X(* (1 \ • • ■ ,- (M) ) (45) 

tlM 
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(cf. (J3J) ) , with the additional restriction that 

X W+1) = 2M + 3 - 2i (46) 
Using the method of derivation of (JsJ) it follows from this that 

P(* (M+1) ) = = n M ( , } (47) 

h m n m =i »^ ! 

which must equal 1 for 

i (M+1) ={l,3,...,2Af + l} (48) 

Consequently 

H M = 2 M(M+1) (49) 

and Caiv = 2^* . Note that 

lim — log A N = 2 lim -|- log Hjvf I 

iV^oc A^ 2 JV^oo A^ 2 ljV/=JV/2-l 

as to be expected from the interpretations of these quantities as entropies for the tiling 
problem. 

There is a second particle system associated with symmetric tilings. This is ob- 
tained by rotating the half Aztec diamond — which has M vertical lines — by 90° 
to obtain a half Aztec diamond positioned with long side horizontal and thus having 
N = 2{M + 1) vertical lines (recall Figure[T|. The first of these is empty of particles 
and last one is full. Ignoring these two lines we have 2M lines where successive lines 
2n— 1 and 2n (n = 1, . . . , M) have n particles. We would like to develop the properties 
of this particle system. 



Analogous to (451, although with Hm substituted by its evaluation (49 1, the joint 



PDF for this weighted particle system is 

P j X *\..., X ™) = ^'V^n (50) 

where x* m the same as x m d5j> , except the first restriction is changed to 1 < x^i < 
x ( p < M + 1. 

Proposition 4.1. Let 

H 2m ^, M = 2 M TT (2i)\ H 2mM = ^) (51) 

i — l 



For p as defined in (50 1, 



/ (2„-r) (2A/), _ A(x<- 2n -V)S(xV n -V) 

P(X '•■■^ j ~ ff 2n _ 1 , M 2»(» (2 "- 1) .-^ (3M - 1) ) 1 j 

Proof. We proceed as in the proof of Proposition |2.1| The 2n — 1 = 1 case is true 



from (50 1. Assume the n = 2m — 1 case is true. Then summing on the (2m — l)-th 
line gives 

.(*») (2m) 
1 «/ (2m-l)\c/.(2in-l)\ 

f/o , .^oM 1 "- 1 ),..,^"^ 1 )) 



(2m-l) (2m) (an. , 



1 det [djA . m (54) 



I7 n , , r o a (^( 2 ™),...^( 2M " 1 )) 1 »>JJ»,i=l,...,m 
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where, with 04 = x^l +2 , 6< = 



This implies 



^ 2m -i,M(2m)!2-(- <2m) ^-^ (2M - 1) ) 



i = 1 
i = 2, 



(55) 



det 



which recalling (51 \ establishes the case n = 2m. Summing now on line 2m gives 

A( JC (2-+D) S ( a; (2-+l)) 



(2m + l) 



, ;a; (2M) ) 



2 m m!# 2 



2 a(*( 3m + 1 ) ) ...,*C3M-l)) 



i,j — 1 , . . . ,m 

(56) 



(57) 



and this recalling (51 1 establishes the case n = 2m + 1. 



□ 



We want to use Proposition |4 . 1 1 to deduce the one line PDFs. For this we again 



a) 

introduce particles y\ representing all lattice sites not occupied by an x particle: 



(n) . . (2M+l-n) 

x K ' U y y ' 

In) _ (2M + l-n) 

x K ' n y ' 



Using the fact that 



we have 



{1,...M + 1} 



(«)! 



5({i,...,Af+i})=n- 



s (as w) = )) n(2y f — ^ - 1) n % M - n+1) n ¥ 

<=i «=i u/j + A -'j-j=i 

where r = |y( 2M_n + 1 '|. It is also not hard to compute that 

(2M + 2)! 1 



(2n) _ U = 

2> 2 2M + 2 (M + 1)! J-jM-n + l^(2(M-n + l)-l) _ lj 



in 

i=l 

Using Proposition |2.2| we have that 

a(x (n) ) =n- Af + a(y (2M_n) ) 



(58) 
(59) 

(60) 
(61) 

(62) 
(63) 



so, changing from x to y in ( 52 I and ( 53 ) respectively we obtain 

jm—l i 



p(y 



(1) „(2m-lh 



p(y (1) ,---,y um; ) 



(2mh 



(M+i-m)! n;L7r( 2 ( M -^+i+ 2 )) ! 

(jVf + 1)! 2 a (y (1, >---^ (2m ~ 2) )+ 2m (- M_m + 1 )+ m 

A(t/ (2m - 1) )g(j/ (2m - 1) ) 

n™x(yf m_1) + M)!(M + 1 - m - x) )! 

nr =1 (2(M-m +J ))! 



2a(3/< 1 ),...,j/< 2m_1 )) + 2m(M-m) 

A(y (2m) )S'(i/ (2m) ) nHi(y' 2m) — |) 



(2m)s, 



(64) 



(65) 
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Using the method of the proof of Proposition |4.1| we compute from these that the 
one-line PDFs are 



(2m-l)\ 



(M + 1 — m)\ (2(M - m + j + 2))! 



(M+l)\ P (2j)! 2 2m ( M ^ m + 1 )+ m 

Afe (2m " 1) ) 2 % (2m " 1) ) 2 
nr=i(yf m_1) + M)!(M + 1 - yf"- 1 *)! 



(2m)", 



(66) 



n 



22m(M-m)-m 



(2(M-m- 



3=1 



1N2 

2 / 



(2m)-., 



(67) 



We remarked above that the one-line PDF ( 24 1 corresponds to a discrete orthogonal 
polynomial unitary ensemble based on a particular Krawtchouk weight. As detailed in 



the Appendix, (64 1 and (66 1 may be regarded as type B versions of the same ensembles. 



We saw in Proposition 



3.1| that the particle system for the Aztec diamond in the 
large N limit relates to the GUE minor process. This is also true of the particle 
system associated with rhombi tiling of the hexagon revised in the Introduction [16] . 
In the case of rhombi tiling of a half hexagon (2a, b, b), cut horizontally along the side 
2a (recall Figure [l]), it is shown in [10] that in the large a limit the particle process 
converges to the eigenvalue process for the minors of anti-symmetric GUE matrices. 
Here we will show that this remains true of the particle system for the half Aztec 
diamond. 

The anti-symmetric GUE is the probability on purely imaginary n x n Hermitian 
matrices with measure proportional to e~ TrX Using the same notation as in (25 1, 
we know from |10) that the joint PDF for the positive eigenvalues of the minor is given 
by 

l/2 n-1 



In- 



c " (=1 

for n even and 



-(*<">> 2 



n (^ n) ) 2 



(4 n) ) 2 



l<j'<fe<n/2 



i=l 



3.(3+1) > x Uh 



1 TT (n) ™') 2 



n 



X 



1=1 

for n odd. Here 



n ((4 n) ) 2 -(4 n) ) 2 )nx(^ 

l<j<fc<(n-l)/2 3=1 



(3+1) > X U)- 



(68) 



, (69) 



Cn 



7r n/4 2-™ /4 , neven 

^(n-l)/4 2 -„(n-l)/4 72 odd 



and it is understood that x 



CO 



0. A straight forward limiting procedure applied to 
(64 1 and (651, according to the strategy of the proof of Proposition 3.1 gives conver- 
gence to these PDFs. 

Proposition 4.2. Let the points z^ 1 ^ := y\^ /y/M be a rescaling of the points y^\ 
where M is the order of the half Aztec diamond as described above, and let z[ = 0. 
Given that the y\ have PDF p as described in (64 1 and 1 65 1, then 



p{y w ,...,y w )^p aGVE {z w ,...,z w ) as N -> oo (70) 



where p a GUE is the PDF for the anti-symmetric GUE minor process specified by 
and Ml 
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Our last task is to compute the limiting support. We proceed using the same 
method as for the full Aztec diamond case. Here however, we will be dealing with 
Boltzmann factors of the form 



n 

l<i<j<N p 



'2 2 1 /3 

X j Xj 



n 



so ( 31 1 becomes, 



V{x) := f 1 P (t) log\x 2 -t 2 \dy. 
Jo 

If we define p(—x) := p(x), then this can be expressed 

V{x) = / p(y)log\x-y\dy 

J —a 



similar to (31 1, although now (32 1 becomes 

p(y)dy = 2N P . 



(71) 



(72) 



(73) 



(74) 



In the limit M —¥ oo, (66 1 and (671 approach a continuum log-gas (71 1 upon the 
substitution 

x^ - - = Mt^ 

where, to leading order, < t < 1. 

In terms of the co-ordinate ti = f- 2 " -1 " 1 , the one body factor in (71 1 reads 



(75) 



-2V0) 



(Af(l+t) + i)!(M(l-t) + i)! 
Let s = n/M. Then A^p = Ms and from ff6| 

'l+t N 



lim — j-j- 4 - = log 



1 -t 



According to (34 1 but taking into account (74 1, a(s) is given by solving 

,- tio g (i±i; 



^ —a 



Va 2 - 



lit = 47TS, 



(cf. (39l) and the integral can be evaluated to give 

^ -2s 



1 - sj\ - a 2 



(76) 



(77) 



(78) 



(79) 



This can be solved immediately for s £ [0, |]. However, for s G (|, 1], like with the 
ear lier case, this equation has no solution for a £ [0, 1]. Using the same logic as leading 
to (41 1, we define a(s) := 1 for s £ (|, 1], so 



a(s) 



Vl-(l-2s 



s £ [0, |] 

sea,i] ' 



giving the same shape as the top half of the full Aztec diamond, as expected. 
In terms of the co-ordinate ti = t[ 2n \ the one body factor in (71 1 reads 



-2V(t) 



Mt 



(M(l + t) + f)!(M(l-i) + A)! 



(80) 



(81) 



This gives the same equation for V'(z) as in ( |77| l, and thus the same result (80 1 for 
a(s), again as expected. 
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A Appendix 



Proposition A.l. Let w(x) be a wetght function with support on successive integers 
a, a+ 1, . . . , b (a > b), and suppose w(x) is even about the midpoint m :— (a + b)/2 so 
that 

w(x — m) — w(—(x — m)), ±(x — m) £ {a, a + 1, . . . , b} (82) 

Let {p n (i)}n=o,i,... be the family of monic orthogonal polynomials p n {x) of degree n, 
with the orthogonality relationship 



}^ W(x)p n {x)p m {x) =(p n ,p n )S 

■m,n 



(83) 



By the property (82), Pj(x — m) — (—l) 3 pj(—(x — m)), so that Pj(x) is even about m 
for j even, and odd about m for j odd. Let w(x) be as in \8£fy and {p n (^)}n=o,i,... be 



as in (83). One has 

b 



" ' w ^ ■ ■■w(x n )(A(x - m) 2 ) 2 = n! J^(p 2j 

x\— a x n —a 
b b n 

H^)---Hxn)(A(x-m) 2 ) 2 l[( 



n-l 

m J_J_tf>2i,P2i) (84) 

3=0 
n-l 

= n\ Yl(p2 3 +i,P2j+il85) 

3 = 



x\=a x n —a i — 1 

Proof. Consider first Q84j> . According to the Vandermonde determinant identity 

A((x - m) 2 ) = det[(a; l - m) 2(j ~ 1) ] ij= i,..., n = det[p 2 (j-i)(a;<)]i,j=i,...,n (86) 

where to obtain the second equality the fact that p2j (x) is even about x — m has been 
used (since A ((a; — m) 2 ) is unchanged (up to sign) by Xi o Xj and by [xi 



— {xi — m) this is referred to as a type B identity). Thus the LHS of (84 1 can be 
rewritten 



^ ••• w(xi)...w(x n )det\p 2( j^ 1 )(xi)]det\p2(j-i)(xi)} 



(87) 



xi-a x n =a 



Since both determinants are antisymmetric in {xj}, while the remaining factors in 
the summand are symmetric, we can replace one of the determinants by n\ times its 
diagonal term Yl7=i P2(i-i)i x i)' We can then perform each sum column-by-column in 
the remaining determinant to reduce (1871) to 



n!det J~] w(x)p 2 (j-i) (x)p 2 (j-i) jx) \ 

L — J i,j=l,...,n 



Now making use of ( 83 1 , (84 1 results 



Regarding ( 84 1, instead of ( 86 1 we use 



A ((a; - m) 2 ) Y\( x i — m) = det[p 2 j-i(a; l )]i J =i, 



which is a consequence of the Vandermonde determinant identity and the fact that 
P2j-i{x) is odd about x = m, and proceed similarly. □ 



15 



The special case p = ~ of the monic Krawtchouk polynomials obey the discrete 
orthogonality equation 

N 1 a! 

(see e.g. [17]). Setting AT = 2M + 1, a; = t + M this reads 
M+i , 

J1 m 2 2M + l (t + M) , (M + ! _ ^P a , 2 M +1 (i+M)p 6 , 2M+ l(t+M) = 22o(2M+1 _ a) ,«5a, i , 

(89) 

If we define a new family of monic polynomials, q n ,N, by q n ,N{x) = p n ,N{N + | + x) 
then we have 

2 2M+i(t + M)!(A./ + l- i )! 9a - 2M+1 (* ~ 2 J 9b ' 2A/+1 (* " 2 j = 2^(2A/ + l-a)!' 

(90) 



The midpoint of the support of the weight in (90 1 is t = 1/2, and furthermore 
the weight is symmetrical about this point. Hence we can apply Proposition |A.1| to 
deduce that 

M+l M+l n 1 , 1 ^2 

E-En (M+I , )!(M+1 _ I / ((- 2 )j 



M+l M+l n _ 1\2 / 1 \ 2 

2 J 



X 1 — 1 X 



= n!C 2n 



where 



&„-, - 22 "''n 2 «(2(A/ 2 -j) + l)! (91) 

*■ - i " , n' i4i ! ^ (92) 

Here the range of summation has been halved by noting Ys x=- m = 2 X^S) 1 wnen the 
summand is even about x = 1/2. Minor manipulation of ( |91[ ) and ( |92| ) reclaims the 
normalizations in (66 1 and (67 1. 
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